Transport equation with nonlocal velocity 
in Wasserstein spaces: 
convergence of numerical schemes 

Benedetto Piccoli* Francesco Rossi^ 
January 31, 2012 

Abstract 

Motivated by pedestrian modelling, we study evolution of measures in the 
Wasserstein space. In particular, we consider the Cauchy problem for a 
transport equation, where the velocity field depends on the measure itself. 
We deal with numerical schemes for this problem and prove convergence of 
a Lagrangian scheme to the solution, when the discretization parameters 
approach zero. We also prove convergence of an Eulerian scheme, under 
more strict hypotheses. Both schemes are discretizations of the push- 
forward formula defined by the transport equation. As a by-product, we 
obtain existence and uniqueness of the solution. 

All the results of convergence are proved with respect to the Wasserstein 
distance. We also show that spaces are not natural for such equations, 
since we lose uniqueness of the solution. 

Keywords: numerical schemes for PDEs, transport equation, evolution of 
measures Wasserstein distance, pedestrian modelling. 
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Recently, when addressing modelling of pedestrian crowd motions, various 
authors used measures to represent relevant quantities, such as pedestrian den- 
sity (see e.g. [3, 4, 6, 7, 9, 10, 12]). The density follows an evolution prescribed 
by a flow map of a velocity function, that is typically composed of two terms: 
a first one called desired velocity, depending on the geometry of the state space 
only, and a second one called interaction velocity, depending on the position of 
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the other pedestrians, thus on the whole measure. The resulting dynamics is a 
transport equation of the type: 



where v = v[fi]. Since each pedestrian interacts with others in a surrounding 
area, the term v is nonlocal, often expressed as an operator depending on a com- 
pactly supported kernel. Moreover, it results very convenient the framework of 
Wasserstein spaces (see [14]), i.e. the space of probability measures endowed 
with the Wasserstein distance. This distance is defined in terms of the solution 
to the optimal transportation problem a la Monge-Kantorovich. Also, numer- 
ical schemes deriving from the push-forward of the measure by the discretized 
flow, proved to be particularly convenient. 

The aim of the present paper is to study convergence of numerical schemes 
for (1) with V = v[fj]. The schemes we study are discretizations of the push- 
forward formula defined by (1). Our main result in this context is that schemes 
in Lagrangian form are the more suitable, converging for any choice of space- 
time discretization. On the other side, Eulerian ones require quite restrictive 
assumptions for convergence. As a by-product, convergence of the schemes is 
used to prove existence and uniqueness of the solution for (1), as proved in [1]. 

All along the article, we use the following notation: we deal with measures 
fj, in A4, the set of probability measures on M" with compact support and 
absolutely continuous with respect to the Lebesgue measure. The first problem 
we address is the existence and uniqueness of a solution for (1) in the case of 
V = v[ij], i.e. w is a vector field depending on the density itself. In this context, 
we always assume the following hypotheses: 



dtn + V • (w/i) = 
M|t=o = Mo, 



(1) 



(H) 



The function 





satisfies 



• V [ij] is uniformly Lipschitz and uniformly bounded, i.e. there exist 
L, M not depending on /x, such that for all ^ G A4, x,y gM", 



\v[ij]{x) -v[ij]{y)\ <L\x-y 



\v[ij] {x)\ < M. 



• V is a, Lipschitz function, i.e. 



there exists K such that 



v[fj]-v [u] II CO < KWp {jj., v) . 
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Under assumption (H), we prove the convergence of a Lagrangian semi- 
discrete in time scheme for (1). We then introduce a complete discrete La- 
grangian scheme, for which we also prove convergence to the solution of (1). 
Given a space discretization parameter Ax and a time discretization parameter 
Af, we prove that the error of the approximation yit^ of the solution on the 
time interval [0, T] satisfies 

Wpip^it),^i{t)) < aAx + bAt 

with a, b given explicitly in Propositions 15 and 16. It is clear that a, b depend 
on the constants L, K, M in (H), the final time T, the metric index p and the 
dimension of the space n. As a consequence, we have weak convergence of fi'^ 
to n for At, Ax — >■ only. 

We then introduce an Eulerian scheme, that has been first proposed in [9]. 
The scheme converges to the solution of (1) under more restrictive hypotheses. 
Indeed, we prove convergence in Wi for At, ^ 0. Wo also prove convergence 

in Wp with p > 1 for At, Ax2^^~^/p^^ — > 0. Also in this case, a precise estimate 
of the error is given in Proposition 18. 

Remark 1 The hypotheses (H) can be relaxed in standard ways. For example, 
the results still hold if we remove the uniform boundedness on R" and ask for 
uniform boundedness only in a point xq. See e.g. [2, Chap. 2]. 

It is interesting to observe that all these estimates are given in terms of the 
Wasserstein distance. This is natural, since in (H) we impose Lipschitzianity 
of V with respect to this distance. For this reason, one could be interested in 
studying the same equation in . This idea is studied in Section 5, where we 
replace Wp with in (H). The surprising result is that we lose uniqueness of 
the solution of (1) under these new hypotheses. For this reason, we don't study 
numerical schemes in this framework. 

The application to pedestrian dynamics also explains the choice of the basic 

assumptions (H), namely that we deal with measures with bounded support. 
The fact that we deal with absolutely continuous measures forbids us to directly 
use the results of this paper for multi-scale models, in which the population is 
modelled by a continuous and discrete part in interaction. Nevertheless, under 
some assumptions (used e.g. in [3]) we can extend the results of this paper to 
that context. This important issue is briefly presented in Section 3.5.2 

The structure of the paper is the following. In Section 1 we study properties 
of the transport equation (1) with v depending on ji. We also recall the deflnition 
of Wasserstein distance and its basic properties. In Section 2 we first study 
properties of the Wasserstein distance under the action of flows. We then study 
the flrst numerical scheme, that is the semi-discrete in time Lagrangian SCHEME 
1. We prove its convergence in C([0, T] , M). As a by-product, we also prove in 
Section 2.3 existence and uniqueness of the solution of (1). 
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The main section of the article is Section 3, in which we introduce the other 
numerical schemes and prove their convergence to the solution under suitable 
conditions. We also discuss in Section 3.5 how to generalize the main results of 
the paper to a finite family of interacting measures (i.e. several populations), 
eventually in a multi-scale modelling. 

In Section 4 we give examples of velocities v = v[ij] that have been in- 
troduced in the context of modelling of pedestrian flows. For the examples 
proposed, we check if (H) are satisfied or not. 

Finally, in Section 5 we discuss the choice of Wasserstein distance or 
distance, both in terms of modelling of pedestrian and in terms of mathematical 
properties. In particular, we show non-uniqueness of the solution of (1) when 
replacing Wp with distance in (H). 



1 Transport equation and Wasserstein distance 

In this section we study the equation (1), in which the velocity v is either given 

by an autonomous vector field v{x) or a time-dependent vector field vt{x). The 
first goal of the section is to properly define the solution of (1) in the context of 
measures. 

Given a measurable map 7 : R" M", one can define the push- forward 
7#/U of a measure /U as follows: 

(7#/x)(i?) = M7"'(J^))- 

The algebra of (7#/i)-measurable sets is thus {E C K" | j~^{E) is /x-measurable} 
The definition of push-forward also implies that, if f = 7#/x, then for all non- 
negative functions cf) it holds 

/ (^07) d/i = / (l)du. (2) 

See a proof in [14, p. 4]. 

Given two positive measures /U, v with the same total mass, one can ask if 
there exists a measurable map 7 such that 7#/U = v. Moreover, among all 
the admissible maps, one can define a cost for 7 of the kind / c{x,'y{x)) dfi 
with c : R" x R" [0, oc) and find the 7 with minimal cost. This is the 
idea of optimal transportation, a problem first proposed by Monge in 1781. A 
complete introduction is given in [14], A particular case is given by the cost 
c{x,y) = \x — y\P with p > 1. This problem induces a definition of distance 
between two measures, called the Wasserstein distance: 

Wpii^, u) = inf I (^J^^ |7(x) - x\P df?j I 7#M = ^ ■ 

This is indeed a distance, see e.g. [14, Ch. 7]. It is important to remark that 
the infimum is always attained, since we deal with measures that are absolutely 
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continuous with respect to the Lebesgue measurc"'^. Another interesting property 
is that the distance can be estimated by dividing the measures in different parts, 
i.e. 

+ A*2, vi + V2) < WP{iii,iyi) + VFP(/i2, 1^2). (3) 

This formula makes sense for positive measures fii, Vi such that = i/j(M"). 

We give two corollaries of this result, that we use in the following. Given /x, u 
sharing a common mass t], we have 

Wj,{ti,v)<Wj,{ix-ii,v-iT). (4) 
Given = Y^f^i and v = Y^'^^i ^^ch that = ^ 0, one has 

Wp{^i,v)< sup Wp{^iuVi)^i{W)-^'P. (5) 

i=l,...,N 

This can be proved as follows. 

N 

W^{li,v) < ^(W;f(/z,,i/,)/Zi(R")-i)/Zi(R")< 

i=l 

/ \p N 

< sup Wp{iii,Ui)fii{R^)-'^p) V/Xi(R") = 



1=1 



sup Wp{iJ,i,iyi)ii^(W') 

\i=l,...,N 

Remark that Wp distances are "ordered", in the sense that 

Pi<P2 implies Wp^{iJ,,i^) < Wp^{fi,u). (6) 

See [14, 7.1.2]. This has a direct consequence in our context. Take a function 

V satisfying (H) for a certain pi with constants L,M,K. Then v satisfies (H) 
for all p > pi and the same constants L,M,K. The converse is not true, since 
we only have that, given X bounded metric space, it holds 



Pi > P2 implies Wp, {p., v) < WP^^p^ (p, u) diarn(Xy- 



P2/P1 



Thus, for p < pi we only have a condition of Holderianity. In particular, the 
strongest condition of the kind (H) is given for p = 1. 

Given a fixed p > 1, we use the following distance in C([0, T] , A4): 

d{n,i/)= sup Wp{nt,iyt)- (7) 

*e[o,T] 

One can also give a Lipschitz family of maps 7* to define the push-forward 
= 7t#M- Under suitable assumptions, one can find a connection between the 
push-forward of a measure and the transport equation. 



^This is one of the main reasons to choose the setting fj, E M. 
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Theorem 2 ([14, Thm 5.34]) Let (7<)jg[ot] ^'^ " locally Lipschitz in time 
family of diffeomorphisms o/ M", with 70 = Id. Let v = v{t,x) be the ve- 
locity field associated with the trajectories of 'jt- Given fio & M., and setting 
IJ't ■= 7t#M0; then fi = /i[o,T] the unique solution of the linear transport 
equation 

dm + V • {vjit) = 

^^\t=o = Mo 

in C ([0,T], At), where M. is endowed with the weak topology. 

The solution of the previous equation is to be intended in the weak sense, i.e. 
for all functions / e ([0, T] x K"), it holds /j„ ^j^jj„ {dtf + Vf -v) diJL = 0. 
We remark that, since we require the continuity of /i w.r.t. i, then it is equivalent 
to state that for a.e. t e [0, T] it holds that, for aU / e C~ (M") 

dt{ I fdfit) = I Vf-vdfif 

See a proof in [14, p. 168]. 

A typical example of an application of Theorem 2 is the case in which v is a. 
given Lipschitz vector field, and 74 = $j is the flow of v. We recall that ^^{xo) 
is the unique solution at time t of 

X — v{x) 
x{0) = xo. 

Then the velocity field associated to 74 is exactly v. One can easily pass to 
time-dependent vector fields vt, assuming that they generate a flow. This is 
verified when v is measurable with respect to time, uniformly Lipschitz in x and 
uniformly bounded. All the results stated above still hold in this case. In this 
context, both approaches (push- forward of measures and transport equations) 
are equivalent. 

Remark 3 Our problem does not flt exactly into the hypotheses of Theorem 2, 
since we endow M with the Wp metric, and the corresponding topology, instead 
of weak topology. 

Nevertheless, observe that we always deal with a compactly supported mea- 
sure Ho, and that the velocity is bounded. Hence, all fXt have compact support, 
contained in supp(^o) + BMt{0). Then change the metric of IR„ outside the 
support of the //,/ to have M" bounded. In this context, Wp metrizes weak con- 
vergence (see [14, Remark 7.13]), thus weak topology coincide with Wp topology 
on M. 

2 Semi-discrete in time Lagrangian scheme 

In this section we introduce a semi-discrete in time Lagrangian scheme to solve 
(1) and we prove that it generates a Cauchy sequence in C{[0,T] ,M.). Since 
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M is complete, then we have a candidate solution for (1); we prove that it is 
indeed a solution. 

We first need some basic estimates for Wp distances under the action of 

flows. Beside being interesting in themselves, these estimates will be useful in 
the following, both for the study of existence of solutions ioi v = v [/x] and the 
convergence of numerical schemes. 

2.1 Wasserstein distance under the action of flows 

In this section, we prove estimates of the distance Wp under the action of flows 
$((•)■ We recall that, given a time t > 0, the flow (•) is a diffeomorphisni of 
M", thus we can see it as a change of coordinates. 

Proposition 4 Let v be an autonomous vector field, Lipschitz with constant L 
and bounded. Let fijV G A4 be two probability measures. Then 

< e^"W^p(M, i^) (9) 

and 

Wp{ii,^:^#i,)<\\v\\cot. (10) 

Proof. Consider a map 7 realizing Wp{ij,, u) = \'y{x) — x\P dii{x) with 7#/U = 
V. Define the map U{y) = ^^{^{^^{y))) and remark that f/#($J'#/i) = 
Thus 

WP{^^,#^,,^^,#u) < [ |C/(t/)-#d($J'#M)(y) = 

= £ ic/($r(^)) - ^n^r I If I di^ui^x^ux)), 

where the last equality is the change of coordinates y ^ x = We 
estimate the Jacobian | ^ |< e^*, see [2]. Hence, 

Gronwall's lemma gives |$j (7(x)) — $t < e^*\j{x) — x\, thus 

WpP($j"#At, < e^*ef^* [ |7(a;) - x\p dn{x) = e^P+^'^^*WP{n, v). 

The second estimate is similar. Take 7(0;) $j (cc) and observe that |7(a;) — 
x\ < ||t;||c7ot. Since 7#/ti = ^ti^^i by definition, then 

WP{fi, < / hix) - x\P dfi{x) < \\v\\P,otP I diiix) = \\v\\%t^. 

Js." Jr" 

□ 

We now turn our attention to an estimate in which the flows are given by 
two distinct vector fields. 
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Proposition 5 Let v, w be two vector fields, both Lipschitz with constant L and 
bounded. Let ijl.,v & M. he two probability measures. Then 

Proof. The proof is similar to the previous one. Consider a map 7 reahzing 
W^{^,u) = Jjj„ |7(a;) — x\P diJ,{x) with 7#/x = v. Define the map U{y) = 

$™(7($'if (y))) and remark that U#{^]'i^fi) = Applying the change of 

coordinates y 1-^ x = ((y), and following the previous proof, we find 

<e^* [ \<^^{j{x))-<^f{x)\Pdii{x) + e^' [ \<^f {x) - '^'^ {x)\p dii{x) . 

^R" 

Gronwall's lemma gives \^'^{j{x)) — < e^*\'y{x) — x\. Define Vrix) = 

\^r{x) - ^t{x)\ and observe that ro(x) = and 

rr{x) < \w{^^{x))-v{^:{x))\< 

< (a;)) - w{^^,{x))\ + \w{^Ux)) - v{K{x))\ < Lr, + \\v - w\\co. 

thus, applying again Gronwall's lemma, we find rt(x) < ll''~^llg° (e-'"* — 1). Hence 
WP{^Ul^,^r*'^) < e(f+i)^*/ \^{x) - x\P d,,ix) + 

We thus have Wp{^f#n,^'^#iy) < e'^^^Wpin,!^) + e^'^^^'-i) _ □ 

Remark 6 These results can be generalized to non-autonomous vector fields 

Wt,vt, if they generate smooth flows. As already recalled, it is verified if they 
are measurable with respect to time and uniformly bounded and Lipschitz with 
respect to space. In this case, in Proposition 5 we have to replace — w\\co 
with sup( \\vt — wtWc"- 

Remark 7 These results can be easily adapted to positive measures that are 
not probability measures but have finite mass, i.e. = C ^ 1. In this case, 

we have 

I^p($J'#M,^r#z^) < e^"VFp(M,t/)+MK")'/^^^^^^^^^— ^11^^- w;||co. (12) 

Remark that the formula is symmetric. Indeed, = j/(R") to have existence 

of the Wasserstein distance W^p(/i, v). 
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2.2 Definition of the semi-discrete in time Lagrangian scheme 

In this section, we precisely define the the semi-discrete in time Lagrangian 
scheme and we prove that it defines a Cauchy sequence in C{[0,T] ,M). The 
idea of the scheme is the following : divide the time interval [0, T] in intervals 
[jAt,{j + l)At]. For each interval, compute the velocity at the initial time 
VjAt = V [MjA*] and use it as a constant on the whole interval, i.e. compute 



SCHEME 1 

Lagrangian, semi-discrete in time, exact velocity 

INITIALIZATION: Fix a time discretization parameter At. 
Take the starting measure fiQ. 

STEP i+1: Given iiiAt, define Vi^t on the interval [iAt, {i + l)At] as 
ViAt •= 'V [MiA*] and compute the corresponding flow ^Jf'"^'. 
For t e [iAt, {i + l)At], define jit ■= ^T^'#HzAf 

STOP: when reaching T. 



Fix a natural number k and divide [0, T] in 2*^ intervals, i.e. choose At = 
Call /i^ the solution of this numerical scheme. We study the convergence of the 
sequence l^fo t] } in C{[0,T] , M.). We prove the following result. 

Proposition 8 Let v satisfy (H), and fiQ be given. Let fJ' = A^fo t] 
structed by Scheme 1 with At = Then, the sequence |MfoT]}^ *^ ^ 
Cauchy sequence in C([0,T] ,M). 

Proof. To simplify the notation, we assume that T = 1. We call := jJ^j 

and := v [m*^] . The corresponding flow is denoted by //'*^ := . Remark 
that we have - v\\\co < KWp{mj,m\). 

We estimate the distance between two successive approximations, that is 
Fix fc e N, i e [0,1] and estimate Wpifi'l , ^t+^). There exists 
j £ {0, . . . , 2*^ - 1} such that t G ^] . Wc estimate Wp{fi'l, (1^+^) in terms 
of Wp{ii''j = Wp{mj,m2j''^), that is the distance at the initial time of 

the interval. We have two cases, as shown in Figure 2.2. 

• t & [^,^]. CaU t' = t-4,. We have /x*' = //;'#m^^ = 
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. t e [|^, Call t' = t - We have = f/4{f^#m% 

= /3^'+^'''+^#(/^^?+^#m^+i). The key difference in these two ex- 
pressions is that, at step k we apply the same wj^, first for ^^tt then for 
t' . At step fc + 1, we first apply v'^'^^ for then v^j'lx for 

We apply (11) and have 

We estimate the first term using the previous case, with t = . For the 
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second term, we have 

- <+i lico < KW,{m'],m',f^,) = KW,{m'], f^i''+'#m'+^) < 

where the last inequality is given by (10). We finally have 

gLt7p(eit' _ 1) / KM\ 
+ -I -[KW,{m'],m'+') + ^j. 

Remark that we have an estimate in terms of Wp{mj j-m^j'^)- Also remark that 
both estimates are increasing as functions of t', thus the maximum is attained 
at the right extreme of the interval. We thus have 



+ -J -[KW,im''j,m',+') + -^ 

= A''Wp{m^j,m^+^) + B\ 

where 

with ifc := ^jfepr- Remark that A'^ and S'^ don't depend on j. 

As a particular case of the last result, we choose t = and have Wp {mj^i , ^2^+2) — 
A''Wp{mj, m2j'^) + B^. Define the sequence for j = 0, . . . , 2^ 

4 := 0, := ^^=4 + 5*= 

and observe that W^p(mf , m^/^) < . We have = B''^^^. Since > 1, 
we also have Wp{mJ^ .m^'^^) < ''^^'^3 = 0, . . . , 2*^. Summing up the results, 

we have Wp(/i( j/ij"*"^) < aSfe for all t e [0, 1], hence < a^k- 

Define the sequence := a!^k. We observe that, if the series ^2^=1 ^ 
convergent, then /x'' is a Cauchy sequence, since d{}j^ ^ jj}) < ^(M"') M''^^) ^ 

X^jC^a-'. We now prove that this series is convergent, using the ratio test 

fc+1 

limfe^.oo 2^ < 1. We have 



B'' (^'=)(2'=)-l A'^+i-l' 
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To compute the limit, wc first notice that Umfe L^+i = 0. Since Lfe+i = Lk/2, 
we have = ^^r^r^- Then 



^fe+l 6^''+^ gPLk+i _ I 2^~^^ 11 1 

lim — — = hm — = — • ? • , , „ = 1 • — • — = — , 

k k e^k eP^f - 1 2*^+2 2 2 4 ' 

A'^ - I ^"^L^ip+i) _|_ eHK[ePLk _ i^(e(*'+^)-'"'= + 1) - 1 

lim . , , ., = lim ■ ^ 



fe ^^=+1-1 fc e2i'fe+i(p+i) + e^±tii£(ePifc+i - l)(e(p+i)^fe+i + 1) - 1 



^2(p+l)Lfc _ 

= lim 



2(P+i)i. _ 1 + K (2 + + 3)L^ + o(Lfc)) (ePi". _ 1) 



T e2(p+i)i'^'+i - 1 + ^ (2 + (p + 3)Lfc+i + o(Lfc+i)) (eP^^+i - 1) 

2{p+l)Lu + o{Lk) + ^^Lk + o{Lk) ii. 
[1X1 — = iini — 

fe 2{p+l)Lk+l+o{Lk+l) + ^^Lk+l+o{Lk+l) fe ife+i 



limlog ((A'^)'') = lim2'=log J^e^^^^P+i) + ^^(e^'^'' - l)(e(P+i)^^ + 1)^ 



= lim2''log(l + 2Lfe(p+ 1) +o(Lfc) +f (l+Lfc+o(Lfc))(pLfc+o(Lfc))(2+(p+l)Lfc+o(Lfc))) 
= lim2''log(l+2Lfc(p+l+f p)+o(Lfc))=limfc2'=(2Lfc(p+l+f p))+2'=o(Lfc) = f (p+l+f p), 

(A^+^f"^' - 1 ep(P+^+TP)+o(i) _ 1 

lim = — TT E^^: = 1- 

fe (Afe)' - 1 ep(^+^+Tf - 1 

We compute limfc = |, thus X^^i Q:*^ < oo, hence the sequence /x*^ is a 
Cauchy sequence. □ 
We also prove the continuous dependence of the approximate solution com- 
puted via the scheme with respect to the initial data. 

Proposition 9 Let v satisfy (H), and iiq^vq he given. Let jj!' and v'^ he con- 
structed using Scheme 1 with At = starting from /xo o,nd uq, respectively. 
Then 

Proof. The proof is similar to the previous one. Fix the following notation 

m'^ := n'^ := v^^^. Then m*^+i = <^^^ "^#m*^ and n^^;^ = 

Take t e [jAi, (j + and call t' = t-jM. Then, using (11) and (H), we 

have Wp{i4,v'^) < e^^*'Wp{m^,n']) + e^*' 1^ (e^*' - l) fWp{m'^,n'^). Since 

the estimate is increasing with respect to t', the maximum is attained at t' = At. 
Using it recursively, we have 
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Since limAt 
2KAt. Thus 




L, then for At sufficiently small, it holds (e 




At 



T 




that is 



T 

A direct computation gives UmAt {e^^*' + 2KM) ^* = e^T+2KT _ 



□ 



2.3 Existence and uniqueness of solution with velocity de- 
pending on the measure 

In this section we prove existence and uniqueness of a solution for (1). The 
proof is similar to the one given in [1]. We first show that, if v depends on the 
measure nt and satisfies (H), there exists a solution of (1). The key idea is to 
use the convergence of the semi-discrete Scheme 1 and to prove that the limit 
is indeed a solution of (1). We then prove uniqueness of such solution. 

We first recall that M. is complete with respect to Wp, see e.g. [13, Thm 
6.18]. As a consequence, C([0,T] , A^) is complete with respect to the distance 
d defined in (7). We now prove that the limit of the semi-discrete Scheme 1 is 
a solution of (1). 

Theorem 10 Let v satisfy (H), and /xq given. Let jjf' = At|) G C([0,r] ,M.) 
be computed using Scheme 1 with At = starting with /xg = Mo- Then the 
limit p, = limfc fj,'' exists and is a solution of (1). 

Proof. First remark that fl exists, since /U*^ ia a Cauchy sequence, by Proposition 
8, and C([0,T] ,M) is a complete metric space. It is clear that JIq = fiQ. To 
simplify the notation, we assume T = 1. 

We now prove that, given a test function / e C^([0, 1] x M"), it holds 



We already know that yJ' computed via Scheme 1 satisfies the following 



for all j = 0, . . . , 2^^ — 1. Remark that in this case v is evaluated &i t = jAt 
only. To prove (14), we prove that 

2'-l /•0+l)Af r /. /■ 




(14) 
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is equal to 0. We have 

2'-l pO + l)At 



(15) < lim^ / / d/Zt(at/ + ^;[/it]-V/) - / rf/z,^ (^J + [m,'^,] • V/) 



j=0 



< 



j=0 



+ l)At 
jAt 



< 



dt 



I 



d{jxt-i4)dtf 



+ 



/ dAtt^ [At] • V/ - / dMt [M^^At] • V/ 



We define Ci := supjQ t^jxR" |V9(/|. If C\ = 0, then f is constant, hence / = 
and (15)= 0. Otherwise, < Ci < oo, then ^ is Lipschitz with respect to 

space with constant at most 1, thus d{flt — l^t)%{' — T^i(A*t)Mt )> using the 
Kantorovich- Rubinstein duaUty formula, see e.g. [14, Ch. 1]. Then recall that 
Wasserstein distances are ordered, in particular Wi(Af,/tit) < Wp{p,t, Ht)- The 
quantity 



/ dptv [fit] • V/ - / dfj,^v [fi'^^t] ■ V/ 



(16) 



can be estimated as follows 



(16) < 



+ 



/ dij.tv [iJ.t] • V/ - / dutv [Mt ] • V/ 

Jr" Jr" 

/ djltv [i^'t] • V/ - / di.'^v [i^l] ■ V/ 

JR" JM" 

/ d/z^ [m,^] • V/ - / di^^v [m,'^,] • V/ 

JR" JR" 

< [ dpt \v [fit] - V [f,^] I |v/| + / d{fit - i^'^)v [^,^] ■ V/ 

JR" JR" 

+ / d„^\v[„'l]-v[f,']^t]\\^f\- 

JR" 

Now recall that [[v [fj] — v [v] [[co < KWp{iJ,, v). Moreover, both v and V/ are 
Lipschitz and bounded, thus v [/z^] • V/ is Lipschitz with a certain constant C2 
not depending on t or k. Thus, again by the Kantorovich-Rubinstein duality 
formula, we have 

(16) < / dtltKWpifit,^i'l)\vf\ + C2W^ifit,^^'l)+ [ d/xti^VFp(Mt,M'At)|v/|. 

Jr" Jr" 
Recalling again that Wasserstein distances are ordered, we have 

(16) < (C3 + C2)VKp(Mt,Mf) + c^Wp{^it^^';^t), 

with C3 := sup[g yjxRn |V/|. Going back to (15), we have 

(15) < lim J2 / dt {Wpifit, ) (Ci +C2 + Cs) + CaWpii^lii'^^t)) ■ 

^* ,=0 "^J^* 
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We estimate the first term by passing to the supremum over t, recalling that 
d(/Lt, u) = supj Wp{iit, i^t)- The last term can be estimated by recalling that 

and applying (10). We then have 

2'-l rU+l)At 

(15) < lim V / dt {d{fl, m'=) (Ci +C2+ C3) + C3M{t - jAt)) . 

Since J-^^^^^\t - jAt) = we have 

(15) < lim (^d(A, m'^) (Ci +C2 + C3) + C^M^^ = 0. 

□ 

We now prove that fxt is Lipschitz with respect to time. 

Proposition 11 Let v satisfy (H), and fj, be a solution of (1). Then /x is 
Lipschitz with respect to time, i.e. 

Wp{nt,lis)<M\t-s\. (17) 

Proof. Since Ht is continuous with respect to time, and v is Lipschitz with 
respect to /i, then v [/it] is continuous with respect to time. Call wt '■= v[^h] and 
observe that it is a non-autonomous vector field, continuous with respect to time. 
Hence the corresponding solution vt of (1) satisfies (10), that is equivalent to 
Wp{vt,Vs) < M\t— s\. By uniqueness of the solution of (1) for non- autonomous 
vector fields, the solution is = /Xj. □ 
We now prove the uniqueness of the solution, as well as the continuous 
dependence on the initial data. 

Theorem 12 Let jj^u two solutions of (1) in C([0,T],A^), with v satisfying 
(H). Then 

WM, v^) < (2e)r*--{2^'8^}l Wp{no, vol (18) 

where \x\ is the ceiling function of x, i.e. the smallest integer not less than x. 
In particular, if /xq = 1^0, then Ht = vt for all t G [0, T]. 

Proof. We first observe that Wp(/if , I't) is Lipschitz with respect to time. Indeed, 
given two times t, s G [0,T], we have 

Wpint, Vt) < Wp{Ht, Ms) + Wpifls, t^a) + Wp{Us, Vt), 

hence \Wp{iJ,t,vt) - Wp{iJ,s,Vs)\ < Wp{iJ,t,lis) + 'Wp{vs,vt) < 2M|t - s|, using 
(17). 
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Wc now prove the continuous dependence on the initial data. Define two 
non-autonomous vector fields f{t,x) := v [nt] {x), g{t,x) := v [ft] {x). Since v is 
Lipschitz with respect to ^ and /U is Lipschitz with respect to time, then both / 
and g are Lipschitz with respect to time and space, thus they generate a flow. 
Applying (11), we have 

pLr /p( Lt _ 1 ^ 

WJ,{^lr,Vr)<e'^^^■'+^^'^Wp{^X<,,V^) + - ^ Slip \\vM-v[Vr']\\co{l^) 

r'e[0,r] 

Define ^(f) := sup.,.£[Q Wp{i^T , t^r) ■ Passing to the supremum over r S [0,t] on 
both sides of (19) and estimating ||w[^t] — w[^'t]||co < KWp{iJ.r,i'T), we have 

hence 

gi,t(p+i)/p 

^ L 

for t such that ^"^ — ^ < 1. We look for a (rough) estimate of this formula. 
We define t' = min {^,87?}- ^^^^ ^^^^^ ^ ~ Ke^' /"^e^' -i) ^ q ^^^^^ 
p > 1, then Lt'/p < W. Since W < ^ < 1, then e"' < 1 + 2Lt'. Thus 
^e"'/y.'_i) ^ K(i+2i:^tO(2LtO < < 1 < 1. We now observe that 

(p{t') < 1^0(0) < 2e(^(0). Using it recursively, wc have (/)(nt') < (26)" 0(0). 

^ 2 

Given a fixed t > 0, then nt' > i for n = \tmax.{2L,8K}~\. Observing that 
Wp{iit,iyt) < (t>{t), we find (18). 

Uniqueness of the solution is a direct application, recalling that Wp{n, v) =Q 
if and only \i ji = v. □ 

We finally prove an estimate about the convergence rate of Scheme 1 to the 
solution. 

Proposition 13 Let v satisfy (H), and given. Let fj, = /U[o,t] be the solution 
of (1), and jj,^ the approximation of fi computed using Scheme 1 with At = 
Let At<pL-^. Then 

d{tx, /) < ^ — >-At. (20) 

Proof. The proof is similar to the proof of Proposition 8. For a time t, let j be 

such that t € [jAt, (j + I) At], and define f/ ^ t — jAt. We estimate Wp{fj4, jJ^) 
with respect to Wp{iJjAt, l^'jAt)- Applying (11), we have 

Wpinu n't) = Wp($f '>#MiAt, ^t'^'hf^'^At) < e(^+^^^*'/^%(M,At, n'jAt) + 

^Lt'/puLt'/p) _i 

+ T 11^ [^^] ~ ^ t^iAt] lie- 

re[jAt,t] 
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We estimate \\v [fir] — v [m^a*] IIc" using (H) and (17), as follows: 

\\V \pr] - V \\C° < K (Wp(/ir, MjAt) + Wp(WAt , M^At)) < 

< Jr(Mr + W-p(M,At,M'At))- (21) 
Coming back to Wp(/Xt, /Uj ) and passing to the supremum in (21), we have 

KMt'e^*'/p(e^*'/p - 1) 

+ 1 

Since the right hand side is increasing as a function of t', one can estimate 
Wp{iJ.t, iJ,^) by putting t' = At. Call 

with i/; := LAt/p. Define the sequence for j := 0, 1 . . . , 2*^ as follows 

Observe that, if t < jAt, then Wp{nt,Ht) < a'j- Since ^'^ > 1, then is 
an increasing sequence, thus for all t E [0,T] it holds Wp(pt, Ht) < a^^, hence 
d{iJL,ii'^) < a^k- We now estimate the value of a^^- E)ue to its definition, we 
have 

^. W^-l ^, i^Me^-(e---l)(-^-"^" + ^^^)^^-^ 

We estimate e^" - 1 < 2Lfc for Lk < 1, i.e. At < pL''^. Then 

ei^(e^- - 1) < e^^^ (^"^IT^) ^ e^^+i)^'' (^^) . 
We also estimate 6^^+'^^^" + > g2Lfc ^ q Summing up, we have 



, KMe^He"--!) (e(^+^)^'- (l + g^))^ -1 

Uok i-il r • 57 ; 

^ L e^^'' - 1 



KMe^" ( e(P+i)^^/P (l ' ^^^^ 
= Ai- 



L(e-f"= + 1) 

For Lfe < 1 we have 1 < 6^*= < e. Moreover, = 2KAt/p, and 2*= log(l + 

2KAt/p) < 2KT/p, hence 

4-= < ^77^ ( e(f+i)ii^/p (1 + 2KAt/pf - 1 ) < At ^ — ^. 

□ 
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3 Lagrangian and Eulerian numerical schemes 



In this section we consider other scliemes for (1), with increasing level of dis- 
cretization. We prove that all these schemes converge to the sohition of (1), 
whose existence and uniqueness has just been proved. SCHEME 2 is only an ini- 
tial discretization in space, while SCHEME 3 is semidiscrete in time with initial 
discretization in space. Scheme 4 is Lagrangian discrete in space and time. We 
conclude definitions of schemes with scheme 5, that is discrete Eulerian. 

Schemes 2 and 3 can be considered as intermediate steps to the definition of 
schemes 4 AND 5, that are the main schemes practically used. Indeed, SCHEMES 
1, 2, 3 are not practically feasible in reality, since they need integration of 
"general" vector fields. Instead, schemes 4 AND 5 can be used since they are 
based on approximations of such vector fields. 

The last part of this section is devoted some simple generalizations of these 
results to equations used in modelling of pedestrian, namely several populations 
in interaction, eventually in a multi-scale approach. 

We start this section by introducing spatial discretization. Fix a space dis- 
cretization parameter Ax > and discretize the space M", dividing it in a grid 
of hypercubes with side Ax. For simplicity, starting from now we will use the 
notation for grids in R^, calling each hypercube simply as "square". This also 
comes from our interest in pedestrian modelling, for which the space is usually 
2-dimensional. We denote the grid with the symbol H. We also use the nota- 
tion XA for the characteristic function of the set A C M". Given a measure n, 
we denote the discretized measure (with an abuse of notation) with the symbol 
H [ij] , computed as follows: given a square □ in the grid H, we define H [/x] on 
this square to be constant and with the same mass as /j. on the square. More 
precisely, we have 

M ^ M 

A*o(n) 



This choice clearly permits to preserve the total measure, i.e. //.(R") = H [fj] (R"). 

We now estimate Wp(/^, H[/u]), using a map 7 that redistributes the mass 
inside each square □ only. Remark that 7 is not optimal in general. We have 
M) < A. hix) - x\P d^i < V^AxP, i.e. 

WpifiMM) < VnAx. (22) 

The same idea can be used to estimate Wp(/U, H [fj]) when is a positive measure 
with mass /x(R") = C ^ 1. We have 

Wp(M, H M) < /i(R")^/^'\/^Aa;. (23) 
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3.1 Scheme with initial discretization in space 

In this section we introduce a scheme with initial discretization in space. The 

scheme is based on a grid H of squares with side Ax. We use the notation /ip-* 

(2) 

to denote the solution of this scheme. We define the initial data /Kq as the 

discretization of the initial data /io, i.e. fx^^ := W\po]- We then compute the 

(2) r (2)1 

approximated solution /j.). ' given by evolution according to v Hi ■ It is the 

exact solution starting from an approximated data. Existence and uniqueness 

(2] (2) 

for (1) give us the good definition of /Xj . Remark that nl is not constant over 
each square □ of the initial grid, and moreover it is not even piecewise constant 
in general. 

We now prove that this scheme is convergent. 

Proposition 14 Letv satisfy (H), and ji^ given. Let jjS-'^^ — m[ot] computed 
via Scheme 2 with parameter Ax, and /i be the exact solution of (1). We have 

< Aa;^^^(2e)r^'"^''^'^•«^>l (24) 
Then, jjL^'^^ converges weakly to fj, in C([0,T] ,M) for Ax 0. 

Proof Using (18), we have that Wp(Aif\/it) < (2e)^*'°'^^^^'^^^Vp(/i^,^\/io), 
thus 

d(^(^),M) < (2e)r^'"-W>Vp(Mp\^o). 

Using (22), we have (24). Since d{iJ,^^\n) converges to for Ax — >■ 0, then /x^^^ 
converges weakly to /i. □ 

3.2 Semi-discrete in time Lagrangian scheme with initial 
discretization in space 

In this section we introduce a third numerical scheme, that is a semi-discrete in 
time Lagrangian scheme with initial discretization in space. 

Fix a space discretization parameter Ax and a time discretization At. We 
perform a first approximation fx^^ ^ of the initial data /xq in space, that is ^q^^ := 
H[/[Xo]- We then compute its evolution using SCHEME 1. The resulting scheme 
is the following. 
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SCHEME 3 

Lagrangian, semi-discrete in time, 
initial discretization in space 

INITIALIZATION: Fix a time discretization parameter At and a 
space discretization parameter Ax. 

,(3) 



Given the starting measure /xq, define := H [^o]- 

,(3) 

and compute the corresponding flow . 

iAf 

STOP: when reaching T. 



STEP i+1: Given fJ^^ti d<^fiiie I'iAi on the interval [iAt, {i + l)Af] as 

ViAt ■= V 



For t G [iAt, {i + l)At], define /^f ^ := 3>^-^*#/xf ^ 



We now prove that also this scheme is convergent. 

Proposition 15 Letv satisfy (H), and jio given. Let ji^^^ = be computed 

via Scheme 3 with parameters Ax and At = and fj, be the exact solution of 
(1). Then 

In particular, /x^^^ converges weakly to n in C([0,T] ,M) for Ax, Ai — ^ 0. 

Proof. Define ly the exact solution of (1) starting from /Jq , that is the approx- 
imated solution of /Lt given by Scheme 2 with parameter h. Thus, /x^^^ is the 
approximated solution of u given by Scheme 1 with parameter At. 

Since d{fi, /i^^^) < d{fi, i')+d{i', fj,^^^), then using estimates (20)-(24), we have 
the result. □ 

3.3 Discrete Lagrangian scheme with velocity of centers 

In this section we introduce a fourth numerical scheme, that can be seen as an 
approximation of the previous one. We indicate with nf the function computed 
via this scheme. The index ^ means "Lagrangian". 

Fix a space discretization parameter Aa; and a time discretization At. We 
perform a first approximation /Iq of the initial data /io in space, that is /i^ := 
H[/xo]. We take a square of the grid and denote it with Dq. We denote the 
center of this square as xa^ and w(no) to indicate the velocity v [/Uq] (xa^), 
i.e. the vector field v [/Iq ] evaluated at the center of the square. This can be 
seen as the approximation of the vector field v [/Xg ] with a vector field that is 
piecewise-constant on the same grid. 
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We compute the evolution Da* of the square Dq as its translation of the 
vector At !>(□□). We define /x^t as the sum of all the translated squares, i.e. 



XDAt 

Ax" 



We repeat the translation □2Af of the square D^t, following v{DAt) v (^nAji 
i.e. the new velocity field evaluated at the center of the squares. We continue 
until reaching T. 



SCHEME 4 
Lagrangian, discrete in space and time, 
velocity at the centers 

INITIALIZATION: Fix a space discretization parameter Aa; and a 

time discretization parameter At. 

Given a starting measure /iqi define as /Xq := H[/io]- 

STEP i+1: Given define for t e [iAt, {i + l)At] 



with Ut ^ DiAt + - iAt) v(niAt)) where v(niAt) := v [/x^At] ip^U^t) 
and XQ.^j is the center of the square DiAt- 

STOP: when reaching T. 



One of the advantages of Scheme 4 is that the evolution of a square is still 
a square, moreover of the same dimension and with sides parallel to axes of M". 
In other words, the measure [i\ is always piecewise-constant. Remark instead 
that the grid is used at the beginning of the algorithm only, and that afterwards 
the squares do not belong to the grid. Moreover, the squares can overlap and 
follow different velocities, thus the evolution cannot be written globally as a 
solution of (1), since v [/i^] (x) is not uniquely determined. 

Another important feature of the scheme is that the value of the function 
inside each square does not change. Nevertheless, there is an interaction among 
all the squares, since the velocity of each square is given by the map v \\x^. One 
can also observe that, if at time ZAt two different squares have the same center, 
then they have the same dynamics from that moment on. This is possible, since 
the evolution of centres is not given by a flow (that would give existence and 
uniqueness), but a discrete-time dynamics. 

We now prove that the scheme is convergent. 



Ax' 



.n 
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Proposition 16 Let v satisfy (H), and /iq be given. Let fi^ = /i^ be com- 
puted using SCHEME 4 with parameters Ax and At ~ Let fi be the exact 

solution of (1). Then d{ii,ii^) < M^^^) + M^) and, for At < ^2^, it 

holds 



d{i,^^\f^^)<2LV^- - 



Ax. 



Then /x^ converges weakly to fj, in C([0, T] , A4) for Ax, At 0. 

Proof. The idea of the proof is to decompose the evolution of measures and 
/zf in tlic evolution of each square of the initial grid. 

Take a square Dq at the starting grid and consider the measure jlo := 
IJ,^a\no- Let t e [jAt, {j + I) At] and call S := jAt, t' = t- S. Consider the 

two following evolutions: the first is the evolution given by /if := #/is- 

The second is the evolution given by := Observe that in both 

formulas, the evolution is computed according to the vector field given by the 



whole measures /Xg'' and fig. As a consequence, we have fi^"^' = J^a i^t ^^'^ 

We denote with xt the evolution of the center of the square Ot according to 
the first evolution, that is Xt = (xs)- The evolution of the center of the 



(3) 



square Dt according to the second evolution is x^ — 



We estimate the distance Wp{jlt,i'f), using the decomposition represented 
in Figure 3.3. 



We have 



r (3) 



We now estimate the four terms. 

We estimate the first term using (12). We estimate \\v 



H'S 



(3) 

Ms 



on the support of jls only. Take a point ys G supp (ps) a-nd estimate \v 

(3) 
Ms 



(^s)llco 

(ys)- 



(3) 

Ms 



(.Ts)|. Since ys G supp (/ts), then it is the evolution of a point yo G 
supp {jls) via Scheme 3, thus 

|yO-+l)At -i0+l)Atl = \yjAt + ^tv flfl^ (yjAt) - (x(^j+i)At + '^tV fxfl^ (XjAt)) I < 
< \yjAt - XjAt \ + ^tL\yjAt - XjAt\, 
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Figure 2: Convergence of Scheme 4: decomposition. 



hence \ys - xs\ = \yjAt - XjAtl < luo - 5o|(l + . Since yo belongs to the 

square Do and xq is its center, we have \yo — xo\ < ^Ax. Recall that v 



is Lipschitz, thus ||w 
up, we have 



(3) 

Ms 



(3) 



(3) 
Ms 



{is)\\c° < L^Ax{l + AtLy. Summing 



l)^Ax{l + AtLy. 



We estimate the second and third terms. Remark that, since a vector field 
evaluated at one point is a vector, then each of these two terms represents the 
Wasserstein distance of the same measure under two different translations. It 
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is easy to prove that, given two vectors a, b, it holds 

< \{x + ta)-{x + tb)\Pd^i 



i/p 



\t\\a- 



This is indeed a particular case of (12). For the second term, we have 

For the third term, recalling the Lipschitzianity oi v, we have 

The fourth term has a particular structure too, since it is the Wasserstein 
distance of two measures under the same translation t' v [f^s] i^s)- Since the 
Wasserstein distance is invariant under translation, we have 

Summing up the four estimates, we have 

Wpifiuv^) < /is(M")'/^e^*'/f (e^*' - l)^Aa;(l + AtLy + (25) 
+t'L\xs - 4|ms(M")'/^ + ^'Wp(4'^Mi)/^s(K")'/^ + Wp{ils,iy^). 

The estimate is increasing as a function of t', thus the supremum is attained 
for t' = At. We now estimate |a;(j+i)At — ^q+i)aJ terms of Wpip^^^^., ^j^^). 
Recall that Xt and x^ are the evolution of the same starting point Xq under two 
different numerical schemes. Thus 



F(j+l)At 



^{j + l)At \ \XjAt + Atv 

< \ijAt-xf^t\ + At{\v 



MjAt 



(XjAt) - V 



(4At)l + 



+ V 



{xf^t)-v[^fAt\ (^iAt)l) < 



< \xjAt - xf^,\ + LAt\xjAt - xf^,\ + KAtWp{^^fl^, fif^,). 

One can observe that |io — Xq\ =0, and also Wp{ii'^\ Hq) = 0, since the two 
schemes have same initial data. Thus \xAt — x^f \ = 0. We thus have 



\XjAt - X 



(1 + Lm) - 1 j£{o,i,...j-i} 
^ ie{o,i,...,j-i} 



We plug this formula in (25). 
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Wc now pass to the global estimate of the distance. Fix a time t and for each 
square □ with /Uo(n) > 0, compute the following quantity: Wp{jlt, 
Now maximise this quantity among all squares and denote it with ^{t). Now 
come back to (25). Wc multiply for /is(M")~-^/P on both sides, then maximise 
on all squares on both sides. It gives 

+ l)At) < e^^*/p(eiAt _ 1) _^. js^tiy + 

+AtK{e^^^' - l)sup,g{o,i,...,,_i}M^p(Mli\,MfAt) + AtKWpij.f , fi^s) + Hj^t). 

Now call fj := sup^^^Q j^^^ (j){t). Since Wp{nl^\ jj^) < 4>{t) due to (5), we have 

fj+i < e^^*/P(e^^* - l)^Ax{l + Atiy + (1 + AtKe^^^')fj. 

We maximize the coefficients (1 + AtLy and e-'^^* with respect to j < 2*^, 
that give e^^ in both cases. We estimate the term under the hypothesis 
LAt < log(2), for which we have e^^* < 1 + 2LAt and e^^*/P < e^^* < 2. 
Since /o = 0, we have 

f . < e-^*/-re-^* - l)^A.e-ii±^^^^^:^l^ < 2LV^Ax'-^^^^—^ 
h^<e [e l)^Axe ^ + ^^^gLT _ i <2WnAx ^ 

Since wc have d{jif-^\ j.i'") < /2fc, we have the estimate. We thus have weak 
convergence of fj,^ to fj,'^^^ for Ax — >■ 0. Using Proposition 15 under the additional 
hypothesis At 0, we have weak convergence of /^^^^ to /i, then of to /U. □ 

Remark 17 The convergence of the Lagrangian schemes is not really surprising 
in the context of pedestrian modelling by using measures. Indeed, we approxi- 
mate a set of discrete agents (the pedestrians) with a continuous measure. This 
is the passage from microscopic to macroscopic model. Then, the scheme is 
convergent for Ax 0, that means that the we go back to a microscopic scale, 
i.e. the pedestrians. 

One can improve this scheme by computing a more precise evolution of each 
square. For example, one can allow deformations of the axes, rotations and 
so on. This idea coincide with the idea of computing the vector field v for 
each square as a certain approximation of the original v. For our scheme, we 
simply evaluate v at the center, i.e. we perform a Taylor expansion of order 0. 
Improvements of this kind certainly result in better convergence rates, but need 
a more complicated implementation. 



3.4 Eulerian scheme 

In this section we present a last scheme to compute numerically the solutions of 
(1). This scheme has been first proposed in [9]. We indicate the solution of the 
scheme as fif. The index ^ means "Eulerian". We call it Eulerian, since we are 
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interested in the evolution of the value iit (P) at a point P not changing in time. 
Instead, the previous schemes where Lagrangian, since we were interested on 
the spatial evolution of a point x, i.e. its trajectory x{t). This was particularly 
clear in the two previous schemes, in which we fixed a starting square Do and 
study its evolution in time 

Fix a space discretization parameter Ax > and a time discretization At. 
We perform a first approximation /i^ of the initial data /io in space, that is 
fiQ := H[/io]- Take a square of the grid □ G H. We denote the center of this 
square as xa and Vt{Oi) to indicate the velocity v [/xf ] (xn), i.e. the vector field 
V [/if] evaluated at the center of the square. 

We compute the evolution of the square □ as its translation of the vector 
Atvo{0), i.e. □ + Atvo{0). We define jlAt as the sum of all the translated 
squares, i.e. yUAt := X]QggMo'(Q) '^°^Aa;"°'°' • ^'^^ moment, this coincide 
with Scheme 4. The difference is that we define /i^^ as the approximation of 
fiAt computing the mean values on the starting grid, i.e. /j.^^ = H[/iAt]- 

We then repeat the same idea "evolution + mean values" starting from 
We contimio until reaching T. 



SCHEME 5 

EULERIAN, DISCRETE IN SPACE AND TIME, 
VELOCITY AT THE CENTERS 

INITIALIZATION: Fix a space discretization parameter Ax and a 
time discretization parameter At. 

Given a starting measure /j-o, define as iIq :— H[/^o]- 
STEP i+1: Define 

XD+Atj)iAt(n) 



M(i+i)At := 2^ MjAtlL-lj- 



Aa;" 



□eB 



with UiAt(n) := [/^^t] (^n) and xu is the center of the square □. 
Define the approximated solution M^+ijAt '■~ ^ [M(i+i)At] • 

STOP: when reaching T. 



We now prove that the scheme is convergent. First observe that is defined 
for times jAt only, thus the distance d{fi,^^) is not defined by equation (7). 
Nevertheless, we redefine it with a slight abuse of notation 

d{fji,iJL^) := sup {Wp(/Xt, /if ) for t € [0, T] such that fxf is defined.} . 

We have the following result. 
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Proposition 18 Let v satisfy (H), and fiQ be given. Let fi^ he computed using 
Scheme 5 with parameters Ax and At = Let ji he the exact solution of (1). 
Then d{iJ,,^^) < d{fi,fj,^) + d{fi^,iJ.^) and we have 



for p > 1: 



d{n\n^) < ^TTTTFTl ^Ax2i^-VP)^. (26) 

For At ^ Q and liniAx.At Ax2^^~^/p'^^ = 0, /u^ converges weakly to fj,; 
• for p= 1: 

For At — ^ and liniAx.At ^ = 0> M'^ converges weakly to ii. 

Proof. We want to estimate Wp{|J,j^^, fif^^). For this reason, we rewrite the 
algorithm in a different way. Assume to have already computed /i^ and /i^ on 
the whole interval [0, T]. Now fix a square Do from the starting grid and follow 
the evolution of the mass inside Do according to the two schemes. Wc denote 
the two evolutions at time jAt with vj" and respectively. Remark that the 
support of lyj" is always a square: we denote such square with Dj" and xf its 
center. The support of Vj" is instead the union of a finite number rUj of squares: 
we denote such squares with , indexed by i = 1 , . . . , rrij . They all belong 



to the grid Effl. We denote with xf^ the center of each square. For technical 



reasons, we always assume that Dj" C U™^■^^Of^, eventually adding squares with 



no mass to the set of squares D^^. 
We have = = ^^Xu,, and 



v. 



J + 1 



with 

rrij \ E ^ E 

^j+i = and := -sit #(1^/ Xni=.). 

One can rewrite the evolutions as follows. Define the vector fields 



v'^ix) := i 



otherwise 



I otherwise. 
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We can rewrite u^^^ = ^^/""^ £^^+1,, = ^"Jf^'^^^^f^i and vf,^^ = H 

Remark that, as already stated, none of the two evolutions is given by a global 
flow, thus we cannot directly apply results of Section 2.1. Nevertheless, take 
7 to be the optimal map realizing Wp{u^ with 'yij^uf = and define 
v^^ := 7#zyjfi. Observe that 

raj 

We estimate the right hand side, that is the Wassertein distance under transla- 
tions of constant vectors v'^(x^) and v^(x^ ^ respectively. Remark that 7 is the 
optimal map realizing Wp{y^^, ^f.i)i by construction. Then it holds 

< [ \y + Atv\xf^,) - (7(2/) + Atv^{xf)))\P duf^M < 

< 2^-1 (^^^ \y - 7(y)r duf^M + At^\v\xf,) - v^ixfruf^M'')) = 

= (l^;(^j;„ ^jf J + AtP\v\xf^,) - v'^{xf)f,.f^Mn) ■ (29) 
Observe that 

\v\xf^,) - v%xf )\ < sup { \v [/if^J (.xf J - V [MiAt] {xf)\ s.t. Df, n Df ^ 0} , 

where we have used that C U^'iDf,-. Remark that Df,- n 9^ when 
jxj' — x^jl < ^/nAx. Hence 

\v\xli) - u"(x^)| < sup {l4M|'At](4)-4A'fAt](^?)l+l''^fAt](xj')-4'"fAt](^j!i) I °j!inn^V0} < 
< if H^p(MiAt, MfAt) + iv^A.T. 

Plug this estimate in (29) and sum over all i in (28). Also recall that, by 
contraction, we have W^{vj',vf) = Yjl=i^p{^j',i^^f,i)- Then compute the 
p-th root, recalling that [a + b)'^/^ < 0}/^ + 6^/^. It gives 

W^{uf+,,^t#uf) < 2^-i/f(w^p(^^,i.f) + AtWp(MtAt,MfAt)^'f(R")'/^ + 

+ AiiV^Axi/f (R")i/f). 

Now observe that 

Wp(i^jVi, i^f+i) < Wp{,^f+„D,+,) + Wp{D,+ui^f+i) < 

< 2i-i/f (Wp{,.^,iyf) + AtKWp{i,f^t,i^f^t)pf(R^y/P + AtLV^Axuf{R^y/p) + 
+x/iAxz/f (R")^/f, (30) 
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where we have estimated Wp{vj+i,vf_^-^) = Wp(yj+i,^{vj+i)) via (23). Mul- 

-i/p 



tiply on the left side for vf_^^ (K") ^'^ and on the right side for vf (K") 
that coincide. Define (/>j := sup m^pi^l ^^f)^i Then (30) reads 



as 

where we have used M^p(/^jAt'A*^t) — ^^^j' proved in (5). Recall that <I)q = Q 
and observe that 2^"^^ + 2^-'^/pMK > 1. Then, for p > 1, we have 

rf(M^,iv^) < sup 0,- < V^Ax{2'-'/PAtL + 1) ^^'"'('^ + trf/^.^? 7 ^ 
^ ^ ~ i=o,i, .,2'= ^-^ - + 2^-^/PAtK - 1 

< VEAx(2^-^/PAtL + 1) — ^-^ < 



< 



< (™^±i)f!IV^A.2»-V..i.. (31) 

- 21-1/P - 1 

This gives estimate (26) and convergence to fj,^ for limAx,At Aa;2(^~^/p)3t = 0. 
The case p = 1 gives 

rf(M^,i.^) < v^Ax(At£ + 1) ^ _ ^ < ^A^(AtL + l)^-^Ax. 

This gives estimate (27) and the convergence to /x^ for limAa;,Af ^ = 0- 

The convergence of /j.^ to n is then given by Proposition 16, since con- 
verges to fl. □ 

Remark 19 The condition given above can be surprising if compared with the 
well-known CFL condition for convergence of schemes for the transport equa- 
tion. We recall that CFL asks for limAx.At ^ > 0) see e.g. [5]. A comparison 
between these two conditions is outside the goals of this paper, nevertheless we 
give two remarks concerning this issue. 

First, it is important to state that Proposition 18 is a partial result, since we 
can't find a counterexample of non-convergence of Scheme 5 under hypothesis 
(H). 

On the other hand, CFL is used for schemes in which the distance between 
measures is the i-'^-distancc. In our setting this is not a natural choice, for 
reasons that are explained in detail in the following Section 5. In other words, 
a scheme that is convergent with respect to one distance can be non-convergent 
with respect to another one. 

Remark 20 It is interesting to observe that, if limAs.At Aa;2*^"'^~^''P^^ = 
is satisfied for a certain pi, then it is satisfied for all p2 < pi- Moreover, 
limAaj.At 'sf = 0) i-e- the condition for p2 = 1 is satisfied too. This is a direct 
consequence of the fact that condition (H) becomes stronger for decreasing p, 
thus the corresponding condition for convergence of the scheme can be weaker. 
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3.5 Generalizations 



In this section we present two simple generalizations of the previous results 
of convergence of numerical schemes. We first focus on the generalization of 
the problem (1) to a problem involving several populations. We then focus 
on problems in which the populations are no more absolutely continuous with 
respect to the Lebesgue measure, but contain Dirac delta too. This second 
generalization is used for multi-scale modelling, see e.g. [3]. 

3.5.1 Generalization to several populations 

In this section we show how to pass from the problem (1) for a population /i 
to the case of a finite number of populations /i^, . . . We explain in the 

following the precise mathematical setting that we are going to study. We 
already observe that this generalization is fundamental for the applications to 
pedestrian modelling, in which one very often has interaction among populations 
with different goals (e.g. pedestrians interacting at a pedestrian crossing). 

We first precisely define the problem. Given iV e N and T > 0, consider 
M := (^1, . . . , /x^) e C([0, T] , each of them satisfying 

dtii' + V-{v' in] M*) = 0, 



In other words, each of the components /i' satisfies its continuity equation as in 
(1), but the vector field [/j,] depends on the global /j,. 

In this context, it is necessary to define distances both for vectors and for 
measures in A^^. Given v,w (i?")^, we define \v — w] := J^^^ — w^\- 
Similarly, given /it, S , we define 

N 
i=l 

Remark that this is one of the possible choices of the distance, corresponding 
to use 1-norm in M.^ . All the following results hold true by changing the norm. 
We also have to change the hypotheses (H) to the following 
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(H-N) 

The function 

v^a^■^ ^ (Ci(K")nL°°(M"))^ 

satisfies 

• V [fj] is uniformly Lipschitz and uniformly bounded, i.e. there exist 
L, M not depending on fi, such that for all /x e A^^, x, y G M", 

\v in] {x) - V in] {y)\ < L\x - y\ \v [n] {x)\ < M. 

• u is a Lipschitz function, i.e. there exists K such that 

\\v[^l]-v[ly]\\co<KWJ,{^l,ly). 

We now review the main results of the paper given for N = 1 under hypoth- 
esis (H). We show that the same results hold for general A'', under hypothesis 
(H — N). We start with the convergence of Scheme 1. The definition of this 
scheme (and of the other ones presented in this paper) is a straightforward 
adaptation of the scheme on each component. 

Proposition 21 Let v satisfy (H-N), and be given. Let /x*' = /^Jj, be 
constructed by Scheme 1 with At = Then, the sequence |a*{ot]} 
Cauchy sequence in C([0,T] ,M^). 

Proof. The proof is completely equivalent to the proof of Proposition 8. Using 
the same notations, one has to estimate \\Vj — V2j~^\\c° < KWp{fj,j , iJi2j^^)- 
Rewrite estimates for each component. Then, summing over all components 
and rearranging the terms, one finds the same recursive formula 

l^p(M^M^') < A'^Wp{m';,m'^+') + B\ 

where and coincide with the ones given in the proof of Proposition 8, 
after replacing K by K N . This gives convergence of the series. □ 
Since C([0, T] ,M.^) is complete, then it exists /x = limfc/x'^. This limit is 
a solution of (32). This can be proved by working on each component as in 

Theorem 10. We have 

Theorem 22 Letv satisfy (H-N), and Hq given. Let = /itfo t] ^ C{[Q,T] ,M^) 
be computed using Scheme 1 with At = starting with /Xg = A*o- Then the 
limit p, = limfe jLt'^ exists and is a solution of (32). Moreover, this solution is 
unique. 
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The uniqueness comes from the same estimate (18), in which constants have 
to be computed more precisely. Similarly, all results of convergence of schemes 
can be adapted to this case, computing the new constants. 

3.5.2 Generalization to multi-scale processes 

In this section we show how to pass from the problem (32) for N populations /i*, 
each of them absolutely continuous, to populations presenting Dirac delta. We 
do not study the general problem in which a Dirac delta can be transformed in 
an absolutely continuous measure. Rather, we consider evolutions in which the 
absolutely continuous and the singular parts cannot melt. In this sense, each 
measure can be always be written as 

/zt = Mr+Mf =7r#/«r+7f#M^, (33) 

where ^1'^ is the absolutely continuous part, and is the sum of the Dirac 
delta. 

This strong condition permits to study this problem as the problem of two 
distinct measures /i"^' and /i'^, each of them satisfying its own transport equation 
(32), with v[ii] — V [/i"'^, /i''] . We recall that the transport equation for a Dirac 
delta Sx is equivalent to the ODE x ~ v{x). Thus one can use Scheme 1 
to find an approximated solution of this problem, then Proposition 21 to check 
convergence, then verify that the limit is indeed a solution. The unique detail to 
verify is that estimates used in the proof of Proposition 21 hold for this special 
case with /it = 7°'^#/io^ + 7^ #Mo- Given this decomposition, one can estimate 
the Wasserstein distance Wp{iJ,,h') by estimating Wp(/i'"^, i^"'^) and Wp{fi'^,u'^). 
Thus, the only result to prove is that estimates given in Section 2.1 also hold 
for measures that are sum of Dirac delta only. A careful look to the proofs show 
that they are verified, mainly using the Gronwall lemma for the ODE x = v{x). 

Finally, it is evident that the same results hold in the case of iV populations, 
each of them satisfying the decomposition rule (33). 

4 Definitions of v for modelling of pedestrian flow 

In this section we present some choices for the definition of v = v [/j] that have 
been proposed for the modelling of pedestrian flows. We briefly describe the 
general idea, referring to [9] for more details. 

Take a pedestrian inside a crowd of people with spatial density /i. He has 
a goal (reach a place, get out of a building, etc.), that defines a vector field 
v'^{x). that is the field of desired velocity. Trajectories to reach the goal are 
thus integral curves of this field. If the pedestrian is alone, his dynamics is only 
determined hy v — v'^. If instead the pedestrian is part of a crowd, it is good for 
him to locally deviate from the trajectory given by f'* to avoid crowded areas 
(that would force him to reduce his velocity). Hence, one can add a velocity 
of interaction that represents the tendency of avoiding crowded areas. It is 
clear that depends on the actual ji, i.e. = [/i] . 
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Thus, a model for the evolution of crowds is given by (1) with 

V = v[ij] = + v'' [ij] . 

We want now study the regularity of v, according to some models that have been 
proposed for v''. As a first observation, the regularity depends on the regularity 
of both terms v''' and v^. We always assume that v"^ is Lipschitz and bounded. 
Under this hypothesis, v satisfies (H) if and only if t;' does. Remark that the 
problem of regularity of v'^ is non-trivial, for example in presence of domains 
with obstacles, see e.g. [8]. 

From now on, we can consider w = w'. Take a pedestrian in x G M." and 
consider how he avoids crowded areas. One can imagine that he checks the 
density of the crowd in a finite area around him and stays away of crowded 
areas. A simple model is the following: take a kernel 77 : K" ^ [0,+oo), 
representing the weight of interaction of the single pedestrian with people in his 
neighborhood. For example, it is reasonable to model that overcrowded places 
close to the pedestrian are "more important" or better detected than far ones. 
Define 

^* lR^yv{x-y)dn{y) 

that is indeed the center of mass of the crowd with respect to the kernel r], and 
define 

V [m] {x) := {x - x*)f r7(x - y) dti{y)^ . (34) 

where / is a non-decreasing weight function. This formula represent the fact 
that V drives away from x* , with an intensity given by /. The meaning of / 
non-descreasing is the following: if the area is crowded, i.e. Jg„ i]{x — y) dfi{y) 
is big, then the pedestrian gets away of the crowd faster than if the area is not 
crowded. Observe that (34) is not defined for 'q{x — y) d^{y) = 0. In this 
case, we simply define v [/i] (x) = 0. This choice does note change the dynamics, 
since an area without mass undergoes no change in time. 

We now study two particular cases of velocity (34). The first is given by 
f{x) = 1, while the second is given by f{x) = x" with a > 1. We prove that in 
the first case v does not satisfy (H), while in the second it does. 

Let f{x) = 1. Given a kernel ry that is not identically 0, we find a family of 
measures /it such that v [/i] is not even a continuous vector field with respect to 
lit in the Wasserstein space. The idea is explained in Figure 4. 

Let R be such that supp (77) C Bfi (Q). Since i] is continuous, we have that the 
set A := {rj > 0} is open. It is always possible to choose r sufficiently small to 
have A\Br{0) nonempty. Since it is open, we can always choose a point x in this 
set, and e > sufficiently small to have B^{x) C A\Br{()). Clearly, if y € B^{x), 
then > r > 0. Finally, define a compact set C of non-zero Lebesgue measure 
outside the ball -Bij(O), and s the maximum distance between elements of Bs{x) 
and C, i.e. s = sup{|a; — y\ s.t. x G B^{x), y e C}. It is clear that s < 00. 
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Figure 3: The convolution kernel does not satisfy (H) for / = 1. 



We are now ready to define the family fit of measures. Define 



A, 



where A is the Lebesgue measure. Observe that v [/io] (0) = 0, since '?(—?;) dfi,Q{y) 
dX{y) = 0. Now observe that, for t > 0, we have /^^ "ni—y)-, dutiv) > 0, 

hence 



> 



inf {|y| s.t. y G B,{S:)} J^^^^^ ri{-y) dX{y) 



> r. 



As a consequence, v [jit] (0) is not continuous with respect to the parameter 
t. We now show that jit is continuous with respect to parameter t in 0, i.e. 

limt Wp(/j,o, M() = 0. Fix a time t and consider the measure vt shared by jiQ and 
fit, that is exactly ut '■= (1 — ^)]^§T ^' Thus, by (4), we have 



'A(C)''''A(i?.(a-))'' 



Take the optimal map 7 between these measures, and observe that ja;— 7(2;)! < s. 
Thus H^p(/Lto, /Xf) < st^/P. Since s and A(C) are finite and do not depend on 
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then Wp{fj.o, i-ii) is continuous at i = 0. Thus, v [fxt] (0) is not continuous with 
respect to the distance Wp{iio, fJ^t)- 

We now study the case f{x) = x" with a > 1. 

Proposition 23 Let v = v[ij] defined by (34), with r] a positive, Lipschitz, 
bounded function with bounded support. Let f{x) = x" with a > 1. Then v 
satisfies (H). 

Proof. Let L be the Lipschitz constant of rj, M = \rj\oo its maximal value, 
R the radius of its bounded support, i.e. supp(r/) C -Bij(O). Call (j){x) := 
/gn vix — y) diJ,{y). We prove that (f) is bounded. Indeed, 



/ 



ri{x - y) dii{y) 



< 



Mdniy) 



= M. 



Remark that this estimate does not depend on the measure fi. 

We now prove that v [/x] is bounded and Lipschitz, uniformly with respect 
to /X. We have 



\v[n] {x)\ 



[ {x - y)r]{x - y) diJL{y) (t>°'~^{x) < ( RM diJ.{y) M""'^ = RM" , 



\v[^i]{x)-v[n]{z)\ 



/ {x - y)v{x - y) di^{y)(l)°' ^{x)+ 



{z - y)T]{z - y) dn{y)(j)°' ^ (z) 



<M"-ii?/ \ri{x-y)-T]{z-y)\dn{y)< 



< M°'-'^RL\x -z\ [ dfi{y) = M°'-'^RL\x - . 



We finally prove that v is Lipschitz with respect to the Wi distance. We have 

a-l 



\v[^i]{x)-v[iy]{x)\ 



/ 7]{x-y)dn{y)] / {x - y)r]{x - y) dn{y)+ 

JR" / 7r» 



^Ir" 



T]{x - y) dv{y) 



{x - y)v{x - y) dv{y) 



< 



/ {x - y)'n{x - y) d{ii - v){y) 

jR" 



since we have that the estimate \ vix — y) dii{y)\ < M does not depend on 
the measure /x. We now prove that f{y) := {x—y)ri{x—y) is a Lipschitz function. 
Indeed, 

\{x - yi)i](x - yi) - {x- y2)'n{x - y2)\ < R\r]{x - yi) - 'q{x - y2)\ < RL\yi - y2\. 

We use this result to estimate J^r,{x — y)ri{x — y)d{ji — i'){y) < RLWidj-jv) 
via the Kantorovich- Rubinstein duality formula. Thus Hf [m] — [j^] ||c° < 
M"-'^RLWi{iJ,,u). □ 
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5 Wasserstein vs distance 

In this section we discuss the choice of the Wasserstein distance Wp as the 
distance in the space M. Indeed, this space is endowed with several other 
distances, see e.g. [11]. A distance which plays an important role in this context 
is the distance, that is 



— h'{x)\dx. 



We will show that the choice of the Wasserstein distance is better than the 

choice of distance, both for modeling and theoretical reasons. 

For the modelling, observe that the Wasserstein metric is more adapted than 
distance to measure if two pedestrian populations are close or far. Indeed, 

take 3 different measures as in figure 4, and call fii the measure centered in Xi, 

for i = 0,1,2. 



I i 



X2 



Figure 4: Definition of /io,/ii,/i2. 



It is easy to prove that 

Wp{^J.i,^I■j) = \xi - Xj\, and - /ijHii =2 for j. 

In particular, the distance between the measures and /^i is identical to the 
distance between fiQ and fj,2- This is not natural with respect to our perception 
of distance between pedestrian crowds, that is better modeled by Wasserstein 
distance. 

Another interesting issue comes from the kind of velocities usually used 
in the context of modelling of pedestrians flows. As showed in Section 4, we 
often deal with velocities defined by convolution with a certain kernel. In this 
case, two measures that are close with respect to the Wp distance but far with 
respect to distance give velocity fields that are close. For example, define 
/Xg :— ■^X(-£.e) s-nd observe that /ig and /i2e are close for small e with respect 
to distance Wp, but not with respect to distance. 
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On the other hand, if one has a velocity v = v[fj] defined by the value of /i in 
a point (provided a good definition of this quantity), i.e. v \p] {x) := f{fj,{x)) for 
a certain /, then the example given above provides a completely different result. 
Choose /(/i(a;)) = We have v [/i^] (0) = 2v [/i2e] (0). Thus velocities are 

far, even if the measures are close with respect to the Wp distance. 

From the mathematical point of view, we show that, in our context, Lips- 
chitzianity with respect to distance is not a good condition for the problem 
(1) with jjQ € ^A. In particular, we show that it does not guarantee uniqueness 
of the solution. Indeed, assume the following hypotheses on v: 



(Hi) 

The function 

v\jj] : 



M Ci(IR")nL° 



satisfies 



• V [fj] is uniformly Lipschitz and uniformly bounded, i.e. there exist 
L, M not depending on /x, such that for all fi G M.,x,y gW^, 

\v [n] {x) - V [ij] {y)\ < L\x - y\ \v [n] {x)\ < M. 

• V is a Li-Lipschitz function, i.e. there exists K such that 

\\vM-v[p] Wco <K\\t,-iy\\Li. 



We now define a function v that satisfies (Hi) and a measure E A4, and 
we show two solutions of (1). The example is based on the classical example of 
non- uniqueness of the solution for the Cauchy problem x = \fx^ a;(0) = 0. 

We first fix the dimension of the space n = 2 and the final time T = 1. We 
define a family of measures i/f S as proposed in Figure 5 (dimensions are not 
respected). 

More precisely, vt is defined as follows. First define a sequence of squares 
Q\, all with sides parallel to axes x and y. Denote with Sj the length of the side 
of the square Q\. All the squares have the upper side on the line y = 1 + t^. 
The square Q\ has the left side on the y-axis. For i > 0, the left side of the 
square Q,\ is contained in the right side of the square Qj-^ Define 



:= ^miXQ'A, 



i=0 

where are positive real numbers to be chosen, and A is the Lebesgue measure. 
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Figure 5: (Hi) does not guarantee uniqueness of the solution. 



We now choose Si := 4~' and nii = ^8*. We first prove that vt G M for all 
t G [0, 1], i.e. it is positive, absolutely continuous with respect to the Lebesgue 
measure, it has bounded support and ft(M") = 1. It is evident that vt is positive 
and absolutely continuous with respect to the Lebesgue measure. We now prove 
that i^t has bounded support. It is easy to prove that the support is contained 
in a rectangle of sides so and J2^o ^i- Since so is finite, we have that i^t has 
bounded support by observing that 

oo oc 14 

i=0 1=0 4 

We finally prove that i't(M") = 1. We have 

oo -. oo -J 

= 5^m,s? = -5^(8/16r = -.— ^ = 1. 

i=Q ^ i=0 ^ 2 

We now define v [i^i] := {0,2t). We first prove that v satisfies (Hi) for i^t 
defined before. We can then extend v to the whole space M so that it satisfies 
(Hi) everywhere. 

For each i>t, v [vt] is a constant vector field, thus v is uniformly Lipschitz 
with L = and uniformly bounded with M = IT = 2. We now prove that v is 
Lipschitz with respect to distance. Let t, s e [0, 1] and assume t> s with no 
loss of generality. Then, for all a; e M" we have 

\v\vt\ {x)-v\v,\ {x)\=Kt-s). 

Take now n such that Sn < — < s„_i. Remark that such a s„ always 
exists, since s„ \ 0. Then, for alH > n the square Q\ is disjoint with respect 
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to Q\. Moreover, all Q\ are always disjoint with respect to Qi for j ^ i. As a 
consequence, 

oo oo ^ 

Pt - z^alUi > Y,'^Vt{Q\) = ^2miS^ = 2m„s2^— — =2-2-", 

i=n i=n 

hence 

||t^[t't]-t^Nllco = 2 (i - s) < 2 - s2 < 27s;;rr = 2 • 2""+^ = 

= 4-2-" < 2||i/t-i/,||ii. 

Thus t; is Lipschitz with respect to distance, with K = 2. 

We have just proved that v satisfies (Hi) for the We now observe that 
we have the two following solutions of (1) with ijlq = vq. The first is the constant 
solution /it = fo for all t € [0, 1], the second is fit = ^t- It is easy to prove that 
the first is a solution, since v [vq] = (0,0), thus constant solutions are solutions. 
For the second, we have to prove that vt satisfies (1) with v [ut] = (0,2<). It is 
equivalent to define a non-autonomous vector field w{t,x) := (0, 2t), consider 
the flow that it generates <I>J"(a;, y) = (.t, y + t^) and observe that ^f^vo = vt- 

It is easy to observe that v does not satisfy (H), for all p > 1. Indeed, 
Wp{vt, Vg) = — , thus it is impossible to find a finite K such that, for all 
t,s, it holds 

\\v [vt] - V [lys] Wco = 2(i - s) < K{t'^ - s^) = KWp{ut, lys)- 

We now extend the definition of v for any ^ G M.. Consider the rectangle 
R := [0,4/3] X [1,2] and define the following functional 

r M ^ [0,+oo) 

Observe that F is well defined and finite, since F(/u) < /u(M^) = 1. We now 
define the function / : [0, 1] — >■ M as follows 

f{F{ut)) := 2t, 

where the Pt are defined above. We have to prove that this function is well 

defined, and defined on the whole domain. It is easy to prove that vt 7^ t^s 
implies F{vt) ^ F{i^s), hence / is well defined. Moreover, F{i't) is clearly a 
continuous function of the parameter t G [0,1], and F(t'o) = 0, F{i'i) = 1. 
Thus, for all ^ G [0, 1], it exists t such that F{i't) = hence / is defined on 
the whole domain. It is also easy to prove that / is strictly increasing. We now 
prove that / is a Lipschitz function. We take < < ^2 < 1 and prove that 
|/(6)-/(a)l <i^l6-ai- We have 

1/(6) 2\t2-h\ 
16-61 |F(z.tJ-F(i.tJ|' 



39 



where ti is the unique value such that Fiyt^) = We first study the particular 



case s 



< tf < t2 < Sn, where Sj 
that all squares Ql with i > n + 1, 



ra+l 



since 1 + 



tk 
> 



= 4 * has been defined above. We have 
k = 1,2, are completely contained in R, 



1. On the other side, each square Ql^ with i < n has 
an intersection with R that is given by the rectangle of horizontal side Si and 
As a consequence, one has 



vertical side 



J2'^tAQlnR)-J2'^tAQk^R) 



i=0 



i=0 



i=0 



i=0 



= 1*2 - ii 1 1*2 + ii I ^ rriiSi > 



i=0 



1 2"+ — 1 

> 1*2 - ill 2,/s;^-J2'^' = Ih - h\ 2-"-i > \t2 - ii|2-"-i2" 



i=Q 



where m, 



8* has been defined above. We thus have 

1/(6) -/(a)i < 2\h-h\ 



16 -a 



\t2-ti\2- 



for satisfying s„_|_i < ii < il < s„ for some n G N. Since the Lipschitz 

constant 4 does not depend on n, we can easily pass to the general case via the 
triangular inequality. One can observe that the Lipschitzianity is also verified in 
0, since / is continuous and the Lipschitz constant 4 for ^ > does not depend 
on ^ itself. 

We now define v [fi] (x) := (0, By construction, we have v [I't] (x) = 

(0, 2f), i.e. this definition of v coincides with the definition given above, that 
gives non- uniqueness. We now prove that this definition of v satisfies (Hi). 
Given /i G A^, then v [/i] is uniformly bounded by 2 and it is uniformly Lipschitz, 
since it is constant. We now prove that v is Lipschitz with respect to the 
distance. Take two distinct measures fi,!/. If they satisfy F{ij,) = F{i/), then 
\\v [fj] — V [v] 11(70 = < ll/z — ^llii. Otherwise, observe that 



IIM"*^!!!/! = / IfJ-ix) — i'{x)\dX{x) > / \iJ,{x) — i'{x)\ d\{x) 



> 



> 



Jr 



(a;)) dX{x) 



\F{f,)-F{u)\. 



Thus 



< 



Hence, v is uniformly Lipschitz if / is a 



Lipschitz function, a fact that we have just proved. 

Remark 24 One can ask if condition (Hi) is sufficient to guarantee at least 
existence of a solution for (1). This is, to our knowledge, an open question. Nev- 
ertheless, the previous discussion is sufficient to show that distance is not 
a good choice for the study of (1) in the case of v depending on the measure itself. 
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